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A brief overview of some recently active topics in perturbative QCD, including: 
string-inspired recursion techniques at tree level; recursion approaches and au- 
tomation of standard techniques for 1-loop calculations; the status of NNLO jet 
calculations; and non-trivial structures that appear in higher-order calculations. 

o 

' 1. Introduction 



As the startup of LHC approaches, much current work in QCD is directed 
towards developing techniques for improving the flexibility and accuracy of 
^H. perturbative calculations. 

<-j ', Flexibility (section 2) is crucial because of the vast range of multi-jet 

final states that will be studied in LHC new-particle searches. At tree level, 
numerical recursion techniques have long been used to build multi-leg am- 
' plitudes from amplitudes with fewer legs — recent developments inspired 

■ by string theory have led to analytically more powerful recursions, giving 

many new compact results for tree-level amplitudes. This new understand- 
ing is also being applied to 1-loop amplitudes, often the missing ingredient 
for quantitatively reliable (NLO) multi-jet predictions. In parallel, more 
traditional 1-loop techniques are being subjected to automation, and here 
too major progress has recently been made. 

Accuracy (section 3), in the sense of NNLO jet calculations, is look- 
ing like it might be within reach in the coming year. This is welcome 
since for over a decade LEP and HERA have been delivering final-state 
measurements with precisions several times better than the NLO theory 
uncertainties, and the latter limit our ability to extract fundamental pa- 
rameters of QCD such as a s and parton distributions. NNLO results also 
provide clues as to the general structures of high orders in QCD. This is 
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both of fundamental interest and potentially useful in predicting large parts 
of yet higher orders. 

2. Multi-jets 
2.1. Tree level 

The multiplicity of Feynman graphs grows factorially with the number of 
legs of a process, hampering the usefulness of traditional techniques for cal- 
culating multi-leg amplitudes. An important discovery of the 1980's was the 
Berends-Giele recursion, 1 allowing amplitudes to be constructed by assem- 
bling smaller, off-shell sub- amplitudes, fig. la. This was suited to recursive 
numerical evaluations and helped prove analytical all-order results, 3 thanks 
to simplifications that occur independently in each sub-amplitude. Berends- 
Giele recursion joins amplitudes via three and four-gluon vertices. Recently 
two new recursions were discovered, CSW 4 and BCF 2 (fig. lb) which join 
amplitudes via a scalar propagator. In the latter case the amplitudes are 
made on-shell using analytic continuation of reference momenta (legs 1,2 in 
figure) . The simpler structure of CSW and BCF recursion makes it easier 
to identify simplifications at each order, leading to many new all-order an- 
alytical results, including applications to processes with heavy quarks and 
electroweak bosons (for a review, see ref. 5). Though originally inspired 
from string theory, the new recursion relations have been proved based on 
general field theoretical arguments 2,6 (exploiting the rationality of tree-level 
amplitudes) and can also be related directly to Feynman diagrams. 7 

The above developments represent an impressive improvement in our 

(a) Berends-Giele recursion 



Figure 1. Graphical comparison of Berends-Giele 1 and BCF 2 recursion relations. 




(b) BCF recursion 
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understanding of field theory. Nevertheless one should bear in mind that 
for practical numerical implementations of tree level calculations, existing 
methods still remain competitive. 8 

2.2. NLO 

While many NLO calculations already exist, 9 there is a recognized phe- 
nomenological need for further multi-leg NLO calculations, in particular 
to simulate backgrounds to new physics signals (for a full discussion of re- 
cent developments, see ref.10). The ingredients in n-leg NLO calculations 
are n + 1 leg tree amplitudes, n-leg 1-loop amplitudes and a procedure for 
combining the pieces. The hardest part is the 1-loop calculation, for which 
several 5-leg results exist and some first 6-leg results are starting to appear. 

The string theory inspired approach to tree-level diagrams does not 
trivially generalise to loop level, in part because of the more complicated 
analytical structure of loop diagrams (cuts as well as poles) . Nevertheless 
the string-inspired approach has led to much new work on loops as well, no- 
tably using the "sewing together" of tree diagrams. This works most easily 
for supersymmetric loop amplitudes, where cancellations between scalar, 
fcrmionic and vector particles in the loop lead to simpler structures in the 
final answer. The plain QCD result is then obtained by combining answers 
with J\f = 4 SUSY, J\f = 1 SUSY and a scalar particle in the loop, the latter 
being the most difficult (split into "cut-constructable" (c, d, e) and purely 
rational (R) pieces). Considerable progress has been made, as illustrated 
in table 1 which shows the contributions to the analytical evaluation of the 
six-gluon 1-loop amplitude, for all independent helicity configurations. For 
MHV configurations (two — helicities, all others +) and split NMHV and 
NNMHV configurations (3 or 4 adjacent — helicities, all others +), general 
multiplicity results have very recently become available. 11 ' 12 

An alternative approach automates traditional methods, i.e. Feynman 
diagram generation, and the recursive reduction of the resulting loop inte- 
grals to a set of known basis integrals. It has the advantage of being easier 
to generalise to processes with external particles other than gluons, but 
suffers from the large number of Feynman diagrams, each term of which is 
broken up into many further terms by the (sometimes numerical) recursion. 
Sometimes the recursion introduces numerical instabilities and alternative 
strategies are then required. 24 ' 25 A notable result with such methods was 
the first full evaluation of the 6-gluon 1-loop amplitude for arbitrary helicity 
configurations, 26 and work is in progress for the 6-quark 1-loop amplitude. 27 
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Full 2^4 NLO jet predictions are however still some way off. 

Related automated methods have been successful also in electroweak 
calculations, with recent full results for e + e~ — > 4 fermions 28 and e+e~ — > 
HHvv, 29 and progress made there will hopefully in part carry over to QCD. 
Also, traditional techniques can simplify considerably 18 when extracting 
just the scalar rational components of the decomposition in table 1 (i.e. 
the parts hardest to obtain in the string-related approaches). 

3. Precision QCD 
3.1. NNLO jets 

Various results exist at NNLO for processes with two QCD partons at Born 
level and one or two non-QCD particles. The current challenge is to address 
processes with three or more QCD legs at Born level, in particular e+e~ — > 
3 jets. All tree- level, 1 and 2-loop amplitudes are known — the difficulty is 
in cancelling divergences between them for a general jet observable. 

Two approaches exist. Subtraction (as at NLO 30 ) identifies a function 
with the same divergences as the real amplitudes, but that is sufficiently 
simple that it can be integrated analytically — one then subtracts the un- 
integrated form from the real amplitudes and adds the integrated form to 
the virtual amplitudes, cancelling all divergences. Finding the subtraction 
functions requires deep understanding of the QCD divergences and ingenu- 
ity so as to make the result integrable. A full scheme at NNLO for processes 
with just final-state particles has been proposed 31 and as a proof of concept 
used to calculate to the cc 3 /^ 3 contribution to the mean thrust in e+e~. 

An alternative approach, sector decomposition, 32 rewrites phase space 
to as to isolate single divergences and then effectively introduces plus- 
prescriptions (as in splitting functions) so as to allow separate extraction of 



Table 1. The analytically derived helicity components of the 1-loop 
6-gluon amplitude (adapted from ref. 13). 
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different powers of the dimension regularisation e. This is less dependent on 
the specific structure of QCD divergences, but becomes more complicated 
as the number of QCD particles increases. It has been successfully used for 
hadron-hadron processes with two Born QCD particles, 33 and for a part of 
the NNLO e + e~ — > 3 jets cross section. 34 

Given the above progress one can perhaps expect first full NNLO predic- 
tions for e+e~ — ► 3jets in the coming year, hopefully with a major impact 
on measurements of the coupling and studies of analytical hadronisation 
models. Extensions to DIS 2 + 1 jet events and hadron-collider dijets will 
probably take somewhat longer. Note that for jets at hadron colliders, an 
issue remains with the experimental jet definitions. Because the standard 
midpoint cone (ILCA 35 ) has the drawback that it can leave large energy de- 
posits unclustered, 36 an extra 'search-cone' step that has been proposed 36 
and used. 37 However this turns out to be infrared (IR) unsafe as the seed 
threshold is taken to zero, 38 compromising theory-data comparisons. A 
positive development is that hadron-collider measurements with the more 
physically motivated (and IR safe) kt algorithm have been shown to be fea- 
sible now by both Tevatron collaborations, 39,40 and the long-standing speed 
issue for the k t algorithm at high- multiplicity has also been resolved. 41 

3.2. Structure of perturbation theory 

Two years have passed since Moch, Vermaseren and Vogt's (MVV) seminal 
calculation of the NNLO splitting functions. 42 With related technology, the 
same authors have obtained the third order coefficient functions, 43 thresh- 
old resummation coefficients, 44 and quark and gluon form factors. 45 These 
results have served as ingredients to calculations of 3-loop Af = 4 SUSY 
splitting functions 46 Drcll-Yan and Higgs threshold resummations, 47 and 
3-loop non-singlet time-like splitting functions. 48 

Various unexpected structures appear in the above results. E.g. writing 

Pij(x) = A + B5(l-x) + C\n(l-x) + 0(l) (1) 
[l x) + 

with A = E™^"( Q! s/47r)", etc., it was noted at NLO 49 that C 2 = A\. At 
NNLO, MVV observed C3 = 2A\Ai. If one postulates splitting functions to 
be universal 50 (identical for time and space-like evolution) when expressed 
for a modified evolution variable z a Q 2 (a = ±1 for the *™° c -like case) a 



Specifically d lnQ2 D = J ^fV{z,Q 2 )D{^ , z a Q 2 ), with D a partem distribution or frag- 
mentation function. 
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and furthermore assumes the universal splitting function to be classical at 
large x (having C = 0), then for normal space-like splitting functions one 
predicts that C = A 2 at all orders, 51 precisely as found at NLO and NNLO. 

The idea of a universal splitting function is given further credibility 
by an analysis 48 which uses the usual 49 analytical continuation x — ► 1/x 
to go from the space-like to the timedike non-singlet (NS) case and finds 
it to be identical to the time-like result found assuming universality with 
the z a Q 2 evolution variable. Note that universality predicts the 3-loop 
P^fg +l — Pns^ 1 difference using only 2-loop information. Given that the 
full 3-loop P^fg +1 and P^fg -1 are themselves also related by x — > 1/x, this 
implies the existence of non-trivial (and yet to be understood) properties 
of the analytic structure of the splitting functions. The universality also 
suggests an explanation for the till-now mysterious absence of 2 and 3-loop 
leading log x terms in the space-like splitting functions, as being closely 
related to exact angular order in fragmentation. 52 Despite these successes 
the universality hypothesis requires further development notably as con- 
cerns the treatment of the singlet sector and the factorisation scheme. 

Other intriguing perturbative structures that have also been found re- 
cently include the following: in Af = 4 SUSY QCD there is increasing 
evidence that n-loop m-leg amplitudes are related to the n th power of the 
1-loop m-leg amplitude 53 (new numerical methods 54 for loop calculations 
providing powerful cross checks); in large-angle soft-gluon resummation for 
2^2 scattering, there is a mysterious symmetry 55 when exchanging the 
kinematic quantity (Ins 2 /ut — 2n) / (ha. u/t) and the number of colours, N c . 

4. Other results 

Owing to limitations of space, many active topics have been omitted. Some 
(small- a; saturation, generalised parton distributions) are reviewed in these 
proceedings. 56 A more extensive bibliography is to be found in ref. 57. For 
others new developments, the reader is referred to the literature, notably 
for 4-loop decoupling relations for a s ; 58 jet definitions that preserve the 
IR safety of flavour; 59 the release of the first C++ ThePEG-based hadron- 
collider Monte Carlo (MC) generator; 60 progress in practical and concep- 
tual aspects of matching MC and NLO; 61 rewcighting to match MC with 
NNLL and NNLO; 62 and soft large-angle resummations, both in terms of 
phenomenology, 63 understanding of treatment of jet-algorithms for non- 
global resummation, 64 two-loop soft colour evolution matrices 65 and other 
recent NNLO resummation results, 66 and an intriguing (but still to be con- 
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firmed) suggestion of a breakdown of coherence at high orders. 67 
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